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ABSTRACT 



This study analyzes the sizes, shapes and correlations of Lya clouds produced 
by a hydrodynamic simulation of a spatially flat CDM universe with a non-zero 
cosmological constant (f2 = 0.4, A = 0.6, er 8 = 0.79), over the redshift range 
2 < z < 4. The Lya clouds range in size from several kiloparsecs to about a 
hundred kiloparsecs in proper units, and they range in shape from roundish, 
high column density regions with Nhi > 10 15 cm -2 to low column density 
sheet-like structures with Nhi < 10 13 cm -2 at z=3. The most common shape 
found in the simulation resembles that of a flattened cigar. The physical size of 
a typical cloud grows with time roughly as (1 + z)~ 3 ^ 2 while its shape hardly 
evolves (except for the most dense regions p cut > 30). Our result indicates that 
any simple model with a population of spheres (or other shapes) of a uniform 
size is oversimplified; if such a model agrees with observational evidence, it is 
probably only by coincidence. We also illustrate why the use of double quasar 
sightlines to set lower limits on cloud sizes is useful only when the perpendicular 
sightline separation is small (Ar < 50/i _1 kpc). Finally, we conjecture that high 
column density Lya clouds (Nhi > 10 15 cm~ 2 ) may be the progenitors of the 
lower redshift faint blue galaxies. This seems plausible because their correlation 
length, number density (extrapolated to lower redshift) and their masses are in 
fair agreement with those observed. 

Subject headings: Cosmology: large-scale structure of Universe - cosmology: 
theory - intergalactic medium - quasars: absorption lines - hydrodynamics 
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1. Introduction 



Acceptable theories of structure formation are those whose parameters have been 
tuned to match the most up-to-date cosmo logical observations. At high redshift, the 
current leading constraint on models is the data from COBE flSmoot et al. 1992[) , which 
fixes the amplitude of the power spectrum on very large scales (~ 1000/i -1 Mpc) to an 
accuracy of about 12%. At low (essentialy zero) redshift, we demand that models fit 
current observations of our local universe, primarily those concerning the distributions of 
galaxies in (x, v) space. These include the abundance of clusters of galaxies, which fixes the 
amplitude of the power spectrum on scales of ~ 8/i Mpc to about 10% accuracy ( Bahcall 
|fc (Jen 1992| ; |Oukbir fc Blanchard 1991 |Bahcall fc (Jen 1993j |White et al. 19934 IViana fc 



iddle 1995| ; |Bond fc Myers 1996| ; [Eke, Cole, fc Frenk 199ft Pen 1996), the power spectrum 



of galaxies, which constrains the shape of the power spectrum on the intermediate-to-large 
scale of ~ 10 - 100/i _1 Mpc (Peacock & Dodds 1994; Feldman, Kaiser & Peacock 1994), 
the ratio of gas to total matter in galaxy clusters, which determines VL b /VL tot ( |White et al.| 



1993b ). In addition, the current measurements of the Hubble constant (Fukugita, Hogan, & 



Peebles 1993; Freedman et al. 1994; Riess, Press, fc Kirshner 1995| ; Hamuy et al. 1995) and 
the age constraint from the oldest globular clusters (Bolte Sz Hogan 1995) limit the range 
for H and the combination of f2 and A . 

This observational suite has been examined by Ostriker & Steinhardt (1995) to 
constrain flat cold dark matter models with a non-zero cosmological constant ( [Peebles 
1984] [Ffstathiou, Bond, fc White 1992| ; |Bahcall fc (Jen 1992| ; |Kofman, Gnedin, fc BahcaTI| 



|1993| |Cen, Gnedin, fc Ostriker 1993| ). The exercise is repeatable for other models (not 
necessarily flat) including the tilted cold dark matter model (|Cen et al. 1992| ; [Liddle,| 
[Lyth, fc Sutherland 1992" ; Lidsey fc Coles 1992| ; |Adams et al. 1993 ; Lucchin, Matarrese, fc 
Mollerach 1993|) , the mixed dark matter model ([Davis, Summers fc Schlegel 1992| ; |Taylor 
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& Rowan-Robinson 1992; |Klypin et al. 1993| ; |Oen fc Ustriker 1994| ; |Ma fc Bertschinger| 



1994| ), the open cold dark matter model (Gott 1982; Bucher, Goldhaber, & Turok 1995) 
and the primeval isocurvature baryon model (Peebles 1987a,b; Cen, Ostriker, & Peebles 
1993). The allowed parameter space for the family of Gaussian cosmological models with 
the tunable parameters (H , Q , A , Rh/c, n i ISO, ADIA) is thus quite limited. In this 
context, H , fl , and A are the Hubble constant, the density parameter, and the value 
of cosmological constant, all in the present epoch. Hh/c is the mass ratio of hot to cold 
matter, n is the asymptotic power spectral index on large scales, and ISO and ADIA denote 
isocurvature and adiabatic models. Although it is possible to further tighten this parameter 
space simply by making more accurate observations of the forementioned quantities, we 
should not overlook any other independent tests which could help to distinguish between 
contending models. 

The Lya forest lines observed in the spectra of high redshift quasars have two unique 
qualities which distinguish them from other observationally accessible phenomena, at 
redshifts between the epoch observed by COBE and our local universe. First, each 
line of sight indiscriminately samples the distribution of neutral hydrogen gas over a 
wide redshift range (z ~ — 5) along random lines of sight (i.e., foreground objects are 
unrelated to the background quasar). Second, each individual spectrum contains a large 
amount of information about low absorption regions (e.g. voids, fluctuating Gunn-Peterson 
absorption) as well as information from high absorption "cloud" lines (number densities 
of clouds at different redshifts, 6-parameters of the individual clouds, correlations of the 
clouds, relationship with other cosmic entities). With a database of thousands of observed 
quasars, the total amount of information available is very large, allowing for very detailed 
statistical studies (e.g., Carswell et al. 1991; Rauch et al. 1993; Petitjean et al. 1993; 
Schneider et al. 1993; Cristiani et al. 1995; Hu et al. 1995; Tytler et al. 1995). These facts, 
taken together, suggest that the Lya forests constitute perhaps the most rich and unbiased 
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sample available for studying the universe at moderate redshift. 

Recent cosmological hydrodynamic simulations by several independent groups have 
consistently shown that Lya clouds are an integral part of the cosmic structure, resulting 
naturally from the gravitational growth and/or collapse of density fluctuations on small-to- 
intermediate scales of ~ 100 kpc to a few Mpc in comoving units (Cen et al. 1994, CMOR 
hereafter; Zhang et al. 1995; Hernquist et al. 1996). Several simple population models have 
been designed to examine the (local) physical as well as (global) statistical properties of Lya 
clouds (Sargent et al. 1980; Ostriker & Ikeuchi 1983; Ikeuchi & Ostriker 1986; Bahcall & 
Spitzer 1969; Arons 1972; Rees 1986; Ikeuchi 1986; Bond, Szalay, & Silk 1988). One of the 
essential simplifications in almost all these models is to assume that individual Lya clouds 
are spherical. While the clouds produced in the new simulations may have some physical 
properties in common with these simple models, a visual inspection of the new simulation 
results reveals that the Lya absorbing structures resulting from small-to-intermediate scale 
structure formation at high redshift are far from spherical. Furthermore, they seem to have 
a wide range of sizes. We presented these results first in CMOR, and then in much more 
detail in Miralda-Escude et al. (1996, MCOR hereafter). Here, we present a quantitative 
study of the topological aspects of the Lya clouds, complementing the topics covered in 
CMOR and MCOR. 

In the redshift range from two to four, Lya clouds exhibit a rich spectrum of structure, 
ranging in shape from semi-spherical to filamentary and even sheet-like. The spherical 
structures often reside at density maxima (with Njji > 10 15 cm -2 ) located in the centers of 
extended structures, while the filaments and sheets tend to have a low Nhi (< 10 13 cm -2 ), 
and form a web-like, interconnecting network covering large portion of the simulation box 
[L = 10/i _1 Mpc comoving). 

We also show that it is not unusual for pairs of quasar sightlines to contain absorption 
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features at close wavelengths. In particular, we show that "common" absorption features 
(coincident lines) in sightline pairs separated by 100/i _1 proper kpc or more are probably 
due to absorption by different clouds, while on separations smaller than 40/i~ 1 kpc, double 
sightlines are mostly likely to actually pierce a common cloud. 

Finally, we show that Lya clouds are spatially significantly clustered, with a 
correlation length of roughly 1 — 2/i~ 1 Mpc comoving (for the high column density clouds of 
Nhi > 10 15 cm~ 2 ) at redshift z ~ 0.5 — 1.0. This suggests the possibility of an intriguing 
connection to faint blue objects (|Koo 1986! [Tyson 1988| ; |Cowie et al. 1988[) . 



This paper is organized as follows. Some brief descriptions of the simulations (for more 
details see MCOR) and our cloud identification method are presented in §2. Results and 
conclusions are given in §3 and §4. 

2. Simulations and Cloud Identification 

2.1. Simulations 

We simulate the formation of Lya clouds in a spatially flat cold dark matter universe 
with a cosmological constant (ACDM), using the following cosmological parameters: 
H = 65km/s/Mpc, tt ,CDM = 0.3645, A = 0.6, tt ofi = 0.0355 (cf. Walker et al. 1991), 
as = 0.79 (the simulations we use in this paper are the same as those used in CMOR 
and MCOR). The primary motivation for choosing this model is that it best fits the 
available observations which were summarized in the introduction. The simulation box 
size is 10/i -1 comoving Mpc per side, and contains N = 288 3 cells and 144 3 dark matter 
particles. The cell size is 35 h~ l comoving kpc, corresponding to a average baryonic cell 
mass of 6.3 x 1O 5 M , with the true spatial and mass resolutions being about 2 and 8 times 



worse than those values, respectively. At z = 3, the Jeans length, Xj = {jic 2 s /Gp 



tot 



,1/2 
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for c s = v rms = lOkm/s, is equal to 400/i -1 kpc in comoving units, or 11 cells. The 
power spectrum transfer function is computed using the method described in Cen, Gnedin, 
& Ostriker (1993). We use a new shock-capturing Total Variation Diminishing (TVD) 
cosmological hydrodynamic code described in Ryu et al. (1993). 

All the atomic processes for a plasma of (H, He) of primeval composition (76%, 24%) 
in mass are included, using the heating, cooling, and ionization terms described in Cen 
(1992). We calculate self-consistently the average background photoionizing radiation field 
as a function of frequency, assuming the radiation field is spatially uniform (i.e., optically 
thin). The evolution of the radiation field is calculated given the average attenuation in 
the simulated box and the emission (both from the gas itself and from the assumed sources 
of ionizing photons). The time-dependent equations for the ionization structure of the gas 
are solved by iteration using an implicit method, to avoid the instabilities that arise in 
solving stiff equations. In general, the abundances of different species are close to ionization 
equilibrium between recombination and photoionization after most of the gas has been 
photoionized. 

We model galaxy formation as in Cen & Ostriker (1992, 1993a,b). The material turning 
into collisionless particles as "galaxies" is assumed to emit ionizing radiation, with two 
types of spectra: one characteristic of star formation regions and the other characteristic 
of quasars, with efficiencies (i.e., the fraction of rest-mass energy converted into radiation) 
of ejjv,* = 5 x 10~ 6 , and euv,Q = 6 x 10~ 6 , respectively. We adopt the emission spectrum 
of massive stars from Scalo (1986) and that of quasars from Edelson and Malkan (1986). 
Details of how we identify galaxy formation and follow the motions of formed galaxies have 
been described in Cen & Ostriker (1993a). Note that in this simulation supernova energy 
feedback into the intergalactic medium from aging massive stars is not included. 



2.2. Lya Cloud Identification 

Perhaps the most critical decision which must be made in a study such as this is how to 
define a Lya "cloud". This is not as simple as it may seem, because the smooth transition 
in density between the global intergalactic medium and the local structures within it can 
blur the distinction between actual clouds and the intercloud medium. However, in order 
to examine the structure of Lya clouds in any quantitative way, it is absolutely necessary 
to adopt some sort of definition of the boundary of a cloud. The simplest approach which 
we believe to be meaningful is to identify clouds as regions with densities above a chosen 
threshold. 

There are two major motivations for defining our clouds in this manner. First of all, this 
method allows us to associate the clouds found at each density cut with a particular column 
density, once we obtain information about the characteristic sizes of the clouds. This is a 
desirable feature because we can only observe Lya forest clouds through absorption lines in 
QSO spectra. Second, density perturbations with amplitudes larger than a certain threshold 
become self-gravitating, bound clouds. For a spherical cloud, the average overdensity at 
which a cloud breaks away from the general Hubble expansion and becomes self-gravitating 
is ~ 5.55. Thus, our cloud defintion also has physical motivation in that structures above 
some high density threshold are gravitationally bound, distinct systems. Having chosen 
a general strategy for defining our clouds, we now describe the details of the grouping 
procedure. 

Cells with a baryonic density below a chosen value (3,10, and 30 will be used) in 
units of the global mean of the baryonic density are cut out of the original density array. 
Then, clouds are defined by grouping the remaining cells using the DENMAX scheme 
( [Bcrtschinger fc Gelb 1991|) as follows. 

First, two sets of baryonic densities are stored. One is the original density output 
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array from the simulation (p org , 288 3 elements), and the second is a smoothed version of 
p org with a Gaussian smoothing window of radius 1.5 cells (p smoo th, 288 3 elements). The 
smoothing is performed to eliminate any small, cell-to-cell density fluctuations, which could 
be physical in the form of small discontinuities or oscillating sound waves, or just numerical 
noises. Such small-scale fluctuations could result in the identification of small unwanted 
clouds with sizes on the order of one cell. At the same time, since the true resolution of 
the simulation is close to 2-3 cells, a Gaussian smoothing window of radius 1.5 cells should 
preserve all information regarding "real" adjacent structures. A larger smoothing window 
could merge real and separate entities which are actually well-resolved in the simulation. 
We have made rather extensive experiments on the smoothing operation and conclude that 
the adopted smoothing window size is appropriate, as will be shown below. 

After the smoothing operation is performed, we cut out all elements of p sm0 oth and p org 
whose p org values are below a chosen value, p cut . We then collect the uncut cells of p org 
(using the DENMAX scheme) by propagating each cell along the gradient of p sm0 oth until it 
reaches a local density maximum. All the cells collected at a particular local maximum are 
grouped into one "cloud" . 

Once the cells have been grouped, we examine the shapes and sizes of the resultant 
"clouds" by approximating them as ellipsoids. A symmetric 3x3 moment tensor is created 
for each cloud of the form T n = Ax 2 ; T 2 2 = Ay 2 ; T 12 = Ax Ay; and so on. This tensor is 
diagonalized, and the square roots of the eigenvalues yield the lengths of the semi-major and 
semi-minor axes, a', b' , c', in descending order of length. For cells where Ax, Ay, or Az = 0, 
we account for the finite size of a single cell by finding /■ , , — = y^ f° r a 1 cen sc i uare an d 
adding that to the tensor component for a particular cloud. This method underestimates 
the actual axial lengths of an ellipsoid, because it only computes moments. We attempt 
to correct this by scaling up the values of a', V ', and d (determined above) by a factor of 
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1/3 



( v 

(a, b, c) = (a', b', c')( %—- ) 1/3 , (1) 

v ' ' ; v ' ' JK 4na'b'c'/3 J ' V ; 

where V is the actual volume of the cloud computed by summing over all the member cells 
of the cloud. 



3. Results 

3.1. A Visual Inspection 

We first present a three dimensional visual description of the structure of the Lya 
clouds produced in the simulation. Figures (la,b,c) display the isodensity contour surfaces 
at z = 3 for p cut = (3, 10, 30), respectively. Here, p cut is the baryonic density in units of its 
global mean. The box size is 10/i -1 comoving Mpc. The most striking visual feature is the 
network of sheets and filaments spanning the box. In Figure la we see that most of the 
structures are connected over the entire box, and it is apparent that most of the covering 
areas [i.e., cross section for quasar lines of sight at low column densities (< 10 13 cm -2 ; see 
below)] are sheet-like objects. Although this result is for the particular ACDM universe 
simulated here, we expect this qualitative picture to be fairly generic for any CDM-like 
model (HDM-like models should show even more prominent sheets, while models like PBI 
are likely to be less coherent). At a higher density of p cut = 10 (Figure lb) most structures 
become filamentary and more isolated. When the isodensity reaches p cut = 30 (Figure lc), 
most of the filaments are replaced by relatively round systems, which are typically isolated 
with separations of about one comoving megaparsec. 

The second noticeable visual feature is the existence of a large, low density region at 
the center of the box, occupying more than half of the simulation volume. Since the initial 
condition of the simulation is not constrained (i.e., it is only a random realization of the 
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cosmological model in such a volume), this void indicates that the simulation box may still 
not be large enough to contain a "fair" volume for the structures under consideration. It 
may be necessary to use a larger simulation box, perhaps 20 — 30/i~ 1 Mpc on a side, in 
order to properly sample the objects in question. Different properties will be affected by 
the simulation volume to varied degrees. We expect that the most significant differences 
will be in quantities like void sizes, or correlations. 

Let us now examine the individual clouds to assess the general accuracy of our 
grouping scheme. Figure 2a shows clouds at p cut = 3 for two randomly selected slices 
of size 5 x 5/i~ 2 Mpc 2 with thickness of I75/i _1 kpc (all lengths are in comoving units) at 
z = 3. The two left panels show the density contours of levels I0^ -1 ^ 2 p cu t,i = 1,2,3..., 
and the two right panels show the identified Lja clouds in the same slices. Two symbols 
are used to show two types of clouds: the filled dots represent "large" clouds with the dot 
size roughly proportional to the actual size of the cloud, and the open circles are for the 
"small" clouds. A cloud is called "small" if its smallest axis does not exceed 3 cells in 
length, and "large" if it does. Figure 2b is similar to Figure 2a, but with p cut = 10, and for 
two different slices. We note that "large" clouds are typically embedded in larger, extended 
density structures while "small" ones are more isolated. A close, one-by-one examination 
of the density maxima in the density contour plots indicates that each identified cloud 
corresponds to a well defined density maximum (Some clouds, especially "small" ones, have 
no corresponding density maxima simply because their densities fall below the contour levels 
when averaged over the slice thickness to be displayed). It appears that our smoothing 
operation and cloud identification scheme indeed yield well defined, distinct clouds with one 
particularly desirable feature: they are neither over-merged in the sense that a single cloud 
contains multiple structures of comparable size (over-smoothed), nor over-separated in the 
sense that artificial clouds are created by small-scale small-amplitude density fluctuations 
(undersmoothed) . 
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When the clouds are classified into "large" and "small" groups, we find that most 
of the clouds are "small": the fraction of "small" clouds at p cut = (3,10,30) at z = 3 
is (62%, 92%, 98%), containing (14%, 49%, 63%) of the total cloud mass. The fact that a 
large fraction of the baryonic mass is in "small" clouds indicates that a higher resolution 
simulation is perhaps needed before we can be absolutely sure that these clouds are resolved 
properly Inclusion of more initial small scale power, which is limited by the Nyquist 
frequency of a simulation, might further increase the fraction of "small" clouds. 

We note in passing that we also experimented with the conventional friends-of-friends 
grouping algorithm, with a linking length of one cell. The resulting clouds were obviously 
not useful for our analyses. For example, at p cu t = 3, most of the volume, as well as mass, 
wound up in a single supercloud because the distinct clouds as shown in Figure 2 were 
linked together by touching boundaries. In other words, the friends-of-friends grouping 
scheme picks out large, inter-connected networks (see Figure la) but not the individual 
structures embedded within them. 

3.2. Quantitative Measures 

We have chosen to focus on three quantitative aspects of the identified Lya clouds: 
their sizes, shapes, and correlations. In addition, double quasar sightline analysis is 
performed. 

To facilitate quantitative discussions, we relate the mass density of a cloud to its 
column density. Assuming that the clouds are in photoionization equilibrium, the column 
density of a cloud can be related to the parallel size and density of the cloud (assuming 
that the density is uniform across the cloud along the line of sight) as 

N H = 8.18 x 10 n ( \-f{J^L-)-\l-)-^— ^_ ( P^ )2( I±f )6 cm - 2 (2) 

v 0.0125/i- 2; v 10- 12; V 10 4; 100kpc(p b ) J v 4 ; ' u 
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where Qf, is the baryonic density parameter, T the temperature in Kelvin, Jhi the hydrogen 
photoionization rate, pb and (pb) the baryonic density of the cloud and the mean baryonic 
density at the redshift in question, L the proper size of the cloud along the line of 
sight, and z the redshift. For a photoionization radiation field with power-law form u~ 1 , 
Jhi = 4.34 x 10~ 12 J LL sec -1 , where Jll is the intensity of the photoionization field at 
the Lyman limit in the usual units (10 _21 erg/cm 2 /hz/sec/sr). Taking fib — 0.0125/i~ 2 , 
Jhi = 7.0 x 10~ 13 sec -1 , T = 2 x 10 4 Kelvin, L = 40 kpc as typical values at z = 3, one 
finds that N HI = (2.5 x 10 12 , 2.7 x 10 13 , 2.5 x 10 14 ) cm" 2 for p cut = (3, 10, 30). Of course, 
exact column densities would depend on exact values of the assumed quantities as well as 
the actual density distribution. We expect that the actual density distribution, which is 
non-uniform, will increase the column densities. For example, redistribution of the same 
mass within L into a (one dimensional) coreless, non-uniform one would result in a column 
density (1 + a) 2 /(l + 2a) times that of the uniform density distribution, where a is the slope 
of the density profile. We simply increase the above estimates by a factor of four to account 
for the gradient in the density distribution of a cloud [note that this factor of four is fairly 
plausible for a reasonable a (~ —1/2) in the one dimensional singular case or a relatively 
lower a value for a non-singular profile. We also note that this factor of four can also be 
achieved if fib is twice the value adopted here, in the light of recent new measurements on 
the deuterium abundance flTytler, Fan, fc Buries 19961 )1, thus obtaining 



N H i = (1.0 x 10 13 ,1.1 x 10 14 ,1.0x 10 15 ) cm- 2 for p mt = (3, 10, 30) at z = 3. (3) 
We will use the above relation for quantitative discussions. 

3.2.1. A Few Simple Global Quantities 

Before engaging in a more detailed quantitative discussion of cloud properties, we 
present a few simple, global quantities. Listed in Table 1 are the numbers of clouds at 
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each redshift and density threshold, the baryonic mass fraction of such clouds, the mean 
inter-cloud separation in proper units, and the mean cloud mass. Two points are interesting. 
First, while the average cloud increases its mass by a factor of about 3 — 4 from z = 4 to 
z — 2, the number of clouds decreases by a factor two. This indicates that the merging of 
old clouds and creation of new clouds are of comparabe importance. Second, the fraction 
of baryonic mass contained in Lya clouds is large. For example, at redshifts (2,3,4), (31%, 
29%, 26%) of baryons are found in regions with densities between p cut = 3 and p cut = 30 
(i.e., N m ~ 10 13 - 10 15 cm- 2 ). 

3.2.2. Sizes ofLja Clouds 
We define the characteristic "size" of a cloud as 

S = (abc) 1/3 , (4) 

where a, b, and c are the semi-axes of the ellipsoid approximating the cloud, measured in 
proper physical units (see §2.2 for definition of a, b, c). For a true ellipsoid cloud, the volume 

is V = f S 3 . 

Figure 3 shows the cumulative distribution of cloud sizes at z = 3 for p cut = (3, 10, 30), 
respectively. The thin curves are weighted by A, where 

A = ab + ac + be, (5) 

and the thick curves are weighted by mass. A is an approximation of the average covering 
area of a cloud on the sky, which is proportional to the probability of the cloud intersecting 
a line of sight. The A- weighted and mass- weighted distributions are very similar, implying 
that A and mass are correlated. Clouds are larger at lower p cut , as expected, and they 
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exhibit a wide range in size, from a few to about a hundred proper kiloparsecs. The median 
size varies from 15/i -1 kpc at p cu t = 30 to 35/i _1 kpc for p cut = 3. Note that S is only the 
cube root of the effective volume of a cloud, so the distribution of cloud sizes would be more 
varied, if we used a (semi-major axis of an ellipsoid) instead of S as a measure of size. 

Figure 4 shows the redshift evolution of the median size of the Lya clouds. Note 
that here the size (vertical axis) is in comoving length units, which might be easier to 
interpret. Our data indicates that Lya clouds grow with time in comoving length units at 
a moderate rate, and that this rate depends weakly on the density (i.e., three curves for 
the three density cases are almost parallel to each other). The fact that denser regions are 
not shrinking, but rather expanding at a rate similar to that of less dense regions implies 
that the merger/accreation rate is higher in the dense regions than in less dense regions. 
Without merger/accretion effects, the high density regions would appear to shrink due to 
gravitational collapse. 

To summarize, we find that Lya clouds exhibit a wide range of sizes from a few proper 
kiloparsecs to about a hundred proper kiloparsecs at z=3. The comoving size of these 
clouds tends to increase slowly with time, implying that the proper physical size increases 
with time more rapidly than (1 + z)~ x . 

3.2.3. Shapes ofLja Clouds 

Figures (5a,b,c) show the distributions of cloud shapes in the c/b - b/a plane for 
clouds at p cut = (3, 10, 30) at z — 3. Each dot represents one cloud and the contours 
indicate the density of clouds (weighted by A) in the plane. The contour levels are 
incremented up linearly from outside to inside. We see that the most common cloud 
at Pcut — 3 and 10 is an ellipsoid with axial ratios of ~ 1 : 2 : 4. However, at p = 30 



-16- 

the situation is interestingly different. Two major concentrations of clouds are seen at 
(b/a, c/b) = (0.35, 0.85) and (0.75, 0.85), indicating the existence of two distinct populations: 
filaments and near-spherical ellipsoids. 

In order to conveniently show the evolution of shapes of the clouds and to relate shapes 
to other quantities, we define a simple "shape" parameter of a cloud as 

b 2 + c 2 
V=-^~. (6) 

For an ideal (spherical, filamentary, disk) cloud, rj = (2,0,1). Some degeneracy exists 
for any value of rj with such a simple definition. For example, (a, b, c) = (1, 1, 0) and 
(a, b, c) = (1, l/v^, 1/V2) both give rj = 1. 

Figure 6 shows rj as a function of size S at z = 3 for p cut = (3, 10, 30). There is a trend, 
albeit with a large scatter, that larger clouds tend to be less spherical. 

Figure 7 represents the A-weighted cumulative distribution of rj for three cases at 
z = 3. We note two points. First, we see that there is negligibly small fraction of clouds 
with the shape parameter 77 greater 1.5, indicating that spherical clouds (rj = 2) contribute 
only a tiny percentage to the total covering area of the clouds at all three density cuts. 
Second, we find a larger fraction of round clouds (large if) and filamentary clouds (77 close 
to zero) at high densities than at a lower densities. This is due to the fact that the large 
sheets and filaments present at lower densities are broken up into both filaments and turn 
into roundish regions at higher densities, respectively. This quantitative result is consistent 
with the visual impression that structures go from sheets to filaments and spheres as the 
density goes from low to high (Figure 1). The fact that there are more filaments than 
spheres at all densities is probably a consequence of nonlinear evolution coupled with 
hydrodynamic (pressure) effects. The situation for dark halos might be different, since dark 
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matter distributions tend to be more unstable in lower dimensional structures. 

Finally, we show in Figure 8 the redshift evolution of the median shape paramter r\ me( i- 
It indicates that at p cu t = 3 the median shape of the clouds do not change with time, while 
at Pcut — 30 the clouds progress toward the shape of spheres (larger rf) in time with a 
dramatic upturn from z = 3 to z = 2, possibly related to the collapse/merging of gas along 
the longest axis. The situation at the intermediate p cu t = 10 is between the above two cases. 

In summary, we have shown that Lya clouds display a variety of shapes ranging 
from semi-spherical clouds to filaments and sheets, with the most common shape being a 
"flattened cigar" with an axial ratio of ~ 1 : 2 : 4 at column density of 10 13 — 10 14 cm -2 , 
and being either a thin cigar with an axial ratios of (1 : 3 : 3.4) or near sphere with an 
axial ratios of (1 : 1.3 : 1.5) at column density of > 10 15 cm -2 . Larger clouds tend to be less 
spherical. The shapes of clouds with densities a few times the mean density of the universe 
evolve very weakly with time, whereas the clouds at higher densities (p cu t > 30) grow more 
spherical with time. 

3.2.4- Double Quasar Line of Sight Analysis 

None of our analyses so far have been directly comparable with observables. Here 
we attempt to project our theoretical models onto the observational plane. CMOR and 
MCOR have examined many of the properties of these clouds in detail, and found that 
the simulation results match observations fairly well. Here, we focus on the observations 
of double quasar sightlines flSmette et al. 1992| ; pinshaw et al. 1994| ; |Dinshaw et al. 1995| ; 



|Smette et al. 1995"t [Fang et al. 1996| Pechtold fc Yee 1996|) . Let us first describe our 



procedure for modeling these observations. 



1. Once three dimensional clouds in the simulation box are identified (see §2.2), each 
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cell is labelled with an integer n, meaning that it belongs to cloud n (if n is zero, it means 
that the cell does not belong to any cloud, i.e., its density is below p cut ). 

2. A random pair of sightlines separated by Ar is selected. Since the simulation box 
has no preferred orientation, we simply choose a direction perpendicular to one of three 
faces of the simulation box as the direction for the double sightlines. 

3. Along each of the two sightlines we identify cells whose labels n (see Step 1) are 
non-zero, and then separate regions of different n into separate clouds, which we call 
Clouds- Along-Line-Of-Sight (CALOS's). Note that each CALOS inherits its 3-d parent 
cloud label n. 

4. We define the location of each CALOS in z-space by finding the p 2 gas weighted center 
of the CALOS cells located along the sightline. 

5. For each CALOS along sightline #1, if a CALOS within a velocity space separation 
D along the sightline #2 is found, the following operation is performed: an integer counter 
N co ,common is incremented if the two CALOS's share the same n, i.e., they belong to the same 
3-d cloud. Otherwise an integer counter N C o, clustering is incremented by 1. If both a common 
and a clustering line are found, we consider it to be a common pairing. The above exercise 
is repeated for each CALOS along sightline #2. We note that the bulk peculiar velocity of 
each CALOS is ignored in this calculation, i.e., the line-of-sight velocity difference between 
a pair of coincident lines simply reflects the line-of-sight real space distance difference. This 
should be a good approximation, because the peculiar velocity gradient is typically smaller 
than the Hubble constant on the scale D (see below) in which we are interested. 

We perform this exercise on 10,000 randomly selected pairs of sightlines for each value 
of pcut = (3, 10, 30) at six different perpendicular separations, Ar = (1,3, 6, 12, 24, 48) cells. 
Figure 9 shows the fraction of lines which are coincident at z — 3, N COytot /N tot , as a function 
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of Ar, where N COttot = N CO)Common + N CO)dustering and N tot is total number of lines. Six 
curves are shown, corresponding to three values of p cut = (3, 10, 30), each at two values of 
D = (50, 150) km/s. We see that the line coincident rate decreases with increasing column 
density. At a perpendicular proper separation of 100/i -1 kpc and p cut = 3 (corresponding 
to clouds with Nhi = 10 13 cm~ 2 ) about 1/2 to 3/4 of the sightlines should have coincident 
absorption features for D = 50km/s to 150km/s. The line coincident rates drop to 
(27%, 38%) and (11%, 17%) for p cut = 10 and p cut = 30 (corresponding to clouds with 
N HI = 10 14 cm~ 2 and N HI = 10 15 cm -2 ), respectively, at D = (50, 150)km/s. 

Figure 10 shows the redshift evolution of N COytot /N tot at Ar = 100/i _1 kpc in proper units, 
for four cases with (N HI ,D) of (3 x 10 13 ,150), (3 x 10 13 ,50), (3 x 10 14 ,150), (3 x 10 14 ,50). 
The indicated neutral hydrogen column densities in Figure 9 are approximately computed 
as follows. Combining Equation (2) for Nhi at z = 3 (with the adopted fiducial values of 
various quantities) with the result that the physical sizes of clouds go roughly as (1 + z) _3//2 
(Figure 4) gives N HI (p cut ,z) = 1.0 x 10 13 (f^) 2 (l + z) 9/2 cm" 2 . We see that, at a fixed 
column density, the rate of coincident lines (at Ar = 100 km/s) increases with redshift at a 
moderate pace, and higher column density clouds and/or smaller D yield weaker evolution, 
assuming that Jhi is constant over the redshift range in question. The assumption of 
constant Jhi is merely for the convenience of illustration. In fact, the self-consistently 
produced photoionization field during the simulation is constant from z = 4 to z = 2 to 
within a factor of two. In other words, the results would not have significantly differed if 
the actual photoionization field had been used. 

Given the three dimensional nature of the simulations, we can distinguish between two 
kinds of coincident lines: one of which (N COyCommon ) occurs when a pair shares the same 
3-d parent cloud, the other of which (N COjdusterin g) occurs when the lines intersect two 
separate 3-d clouds. The former set of coincident lines is due to the extended size of a 
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single cloud, while the latter set of coincident lines is due to the clustering of several clouds. 
Figure 11 shows N COjCommon /N co>tot at z — 3. Six curves are shown for p cut = (3,10,30) 
with D = 50 km/s and 150 km/s. We see that at Ar m 30 — 60/i _1 kpc about half of the 
coincident line pairs share the same clouds and the other half pierce separate clouds. At 
Ar = 100/i _1 kpc only 10-20% of the total coincident line pairs share the same clouds, i.e., 
clustering of separate clouds dominates the coincident events at this and larger separations. 

The observational signature of clustering dominated coincident lines is that the 
difference in line of sight velocity between the two absorption lines should be a weak 
function of the separation perpendicular to the line of sight. Figure 12 shows rms velocity 
difference of coincident line pairs along the line of sight, R\\ rms , as a function of the 
perpendicular separation, Ar, for D = 50 km/s and 150km/s at p cut = (3, 10,30). The two 
long-dashed horizontal lines show the rms difference for the two D values, if clouds are small 
(compared to perpendicular separation) and randomly distributed (i.e., no clustering). We 
see that R\\ trms depends rather weakly on Ar when Ar > 200/i~ 1 kpc, and it approaches the 
asymptotic limit of a random distribution of clouds at Ar f» 400 — 500/i -1 kpc. We predict 
that for quasar double sightlines at large separations (Ar > 500/i -1 kpc), the rms velocity 
difference between absorption features common to both sightlines, R\\ trms , should depend 
on the range being searched in z-space, D, as 

R llrms = 0.58D. (7) 

Figures (13a,b) show the redshift evolution of N COiCommon / N COjtot for Nhi = 3 x 10 13 cm~ 2 
and Nhi = 3 x 10 14 cm -2 , respectively. Four cases are shown in each figure with two 
separations, Ar = 50/i _1 kpc and Ar = 100/i _1 kpc both in proper units, and with 
D = 50km/s and 150km/s. We see that the fraction of coincident lines with a separation 
of 50 — 100/i -1 proper kpc due to common clouds decreases with increasing redshift at fixed 
column density, presumably because clouds are smaller and more crowded at high redshift 
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than at low redshift. We remind the reader that the total line coincident rate increases with 
redshift (see Figure 10), indicates that a higher cloud number density (crowding of clouds) 
at high redshift dominates the line coincident events. This is consistent with results shown 
in Figure 11. 

The combined results of our quasar double sightline analysis suggest a few basic 
trends which should be kept in mind when analyzing QSO spectra. The most important 
of these is that line coincident events can be caused by two different phenomena. At 
small perpendicular separations (Ar < 40/i _1 kpc proper) a pair of coincident lines is most 
likely to pierce a common cloud, as is usually assumed when analyzing such observations 
to infer the actual cloud size. However, at larger Ar, a pair of coincident lines is more 
likely to penetrate two different clouds which are spatially clustered. At Ar = 100/i -1 
kpc, only 10-20% of the total coincident line pairs share the same clouds. And at very 
large separations of Ar > 500/i _1 kpc, the coincident line events can be explained as the 
random intersection of two unrelated clouds, whose rms difference R\\ yrms should relate to 
the velocity interval D as indicated in Equation (7). 

3.2.5. Correlations ofLja Clouds 

It is tempting to make a connection between Lya clouds and dwarf galaxies and/or 
moderate redshift faint blue galaxies. We approach this problem by examining clustering 
properties of the Lya clouds. 

A few cautionary words about the limitations of the simulation are in order here. The 
simulation box size (L = lO/i^Mpc comoving) places severe limits on our ability to study 
the clustering properties of clouds on large scales with a high degree of accuracy, because 
it sets an upper limit on the scale of the input power spectrum. Waves longer than the 
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simulation box size would have made a considerable contribution to clustering on larger 
scales, since the density fluctuations on scales comparable to the simulation box length have 
started to approach nonlinearity even at redshift z = 2 [a ~ 0.4 at z = 2 on the box scale). 
The situation becomes even more serious at lower redshift. 

Nevertheless, we may gain some insight into the distribution of mass on smaller scales, 
and say somthing useful about the more local clustering properties of clouds. The three 
dimensional correlation of the Lya clouds at z = 3 is shown in Figure 14 for p cut = (10, 30). 
Higher column density clouds seem to be more strongly clustered than lower column density 
ones. Since most of the faint blue galaxies are thought to be in the redshift z ~ 0.5 — 1.0, 
we show in Figure 15 the evolution of correlation length r (solid curve). This is defined as 
the length where the correlation is unity and is plotted for the redshift range z = 2 — 4, 
where we have the most confidence in the accuracy of the simulation. All lengths are shown 
in comoving units. 

In order to estimate the effect of missing power at box-sized scales, we computed 
the correlation function of the clouds in eight subboxes. These subboxes were created by 
dividing the original simulation box into eight equal cubes and pretending that each subbox 
also has periodic boundary conditions. The dotted line in Figure 15 shows the median value 
of the correlation length for the subboxes, and the vertical bars indicate the full width of 
the distribution at each redshift (i.e., highest and lowest values of r ). We see that the 
median correlation length drops about 30-40% for a box of size L = 5/i -1 Mpc. Our best 
estimate is that the true correlation length at each redshift would be larger than values 
indicated by the solid curves by perhaps 50%, were the simulation box sufficiently large. 

Simple extrapolation of the solid curve in Figure 15 to lower redshifts is unlikely to 
give us accurate values of r at z ~ 0.5 — 1.0. However, it is reasonable to expect the value 
of r at these redshifts to be significantly larger than 1.5/i _1 Mpc comoving, when one takes 
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into account the effect of the missing longer waves. This correlation length is interestingly 
close to the observed correlation length of the faint blue galaxies QEfstathiou et al. 1991| ; 



Neuschaefer, Windhorst, fc Dressier 199 It |Couch, Jurcevic, fc Boyle 1993| ; |Roche et al 



1993| ; [Infante fc Pritchet 1995| ; [Brainerd, Smail, fc Mould 1995|) , although an accurate 



conversion from observed angular correlation function w(8) to the 3-d correlation function 
as computed here requires detailed knowledge of the redshift distribution of the observed 
faint blue objects, which is currently unavailable. Furthermore, the average mass of these 
objects is close to 10 9 — 10 10 M (see Table 1), which is in accord with what the mass of the 
faint blue galaxies is thought to be or required [cf. , e.g., Babul & Rees' (1992) starburst 
dwarf galaxy model]. Finally, the number density of the clouds (with a mean separation of 
~ 2/i~ 1 Mpc at z=2 for p cut = 30) and its increasing trend with redshift indicate that it is 
probably quite close to what is thought to be for observed faint blue objects. These three 
considerations lend tentative support to the conjecture that the high redshift Lya clouds 
(Nhi > 10 15 cm~ 2 ) are the progenitors of the faint blue galaxies. 

4. Conclusions 

We have presented a quantitative study of the sizes, shapes, and correlations of Lya 
clouds in a spatially flat cold dark matter universe with a cosmological constant, utilizing 
state-of-art cosmological hydrodynamic simulations including detailed atomic physics for 
a plasma of primordial composition. Keeping in mind the unavoidable limitations of 
such numerical experiments due to limited boxsize and limited resolution as well as some 
approximate treatment of physical processes, a few findings are probably fairly reliable. 

Structures formed at high redshift z ~ 2 — 4 due to gravitational growth/collapse 
of cosmic density perturbations at small-to-intermediate scales (lOOkpc to a few Mpc 
comoving) are responsible for the observed Lya forest. These Lya clouds cannot be 
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described by a simple model which characterizes their sizes and shapes. Their sizes vary in 
a wide range from a few kiloparsecs to about one hundred kiloparsecs in proper units with 
the median size being ~ 15 — 35/i _1 kpc at z — 3, and their shapes vary from nearly-spherical 
ellipsoids to filaments and pancakes. A typical Lya cloud resembles a "flattened cigar," 
with an axial ratio of ~ 1 : 2 : 4 at column density of 10 13 — 10 14 cm -2 , and is either a thin 
cigar with an axial ratios of (1 : 3 : 3.4) or near sphere with an axial ratios of (1 : 1.3 : 1.5) 
at column density of > 10 15 cm -2 . Larger clouds are, on average, less spherical than smaller 
ones. The physical size of a typical cloud grows with time roughly as (1 + z)~ 3 ' 2 while its 
shape hardly evolves (except for the most dense regions p cu t > 30 where tend to be more 
spherical with time). 

Analysis of simulated quasar double sightlines indicates that coincident absorption 
features observed in the two spectra can have two different causes. At small perpendicular 
sightline separations (Ar < 40/i~ 1 kpc proper), a pair of coincident lines most likely 
represents absorption from a common cloud. This is usually assumed when analyzing such 
observations to infer the actual cloud size. However, at larger Ar, a pair of coincident 
lines most likely samples two separate clouds which belong to the same cloud cluster. 
In other words, clustering of separate clouds dominates the coincident events at larger 
double sightline separations. At Ar = 100/i -1 kpc, 80-90% of coincident line pairs can 
be explained in this fashion. At very large separations, Ar > 500/i -1 kpc, the coincident 
line events are entirely due to random intersections of two unrelated, uncorrelated clouds, 
whose rms velocity difference R\\, rms should relate to the velocity interval being examined, 
D, as R\\ rms = 0.58-D. Analyzing observed coincident absorption lines to infer the actual 
sizes of Lya clouds by assuming a population of clouds with simplified geometry (e.g. 
spheres, disks) is reasonably accurate only when the double sightline separations are small 
(Ar < 50/i _1 kpc). For larger Ar, this exercise is not very meaningful. The inferred sizes 
are not related to true sizes of the clouds and are grossly inflated. 
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We attempt to make a connection between the high column density Lya clouds 
{Nhi > 10 15 ) and the faint blue galaxies. The supportive evidence for such a conjecture 
is threefold. First, when extrapolated to z ~ 0.5 — 1.0, the correlation of the high column 
density Lya clouds (> 1.5/z _1 Mpc comoving) is close to that of the observed faint blue 
galaxies QEfstathiou et al. 1991! [Neuschaefer, Windhorst, fc Dressier 1991| ; |Couch, JurcevicJ 



k Boyle 1993| ; |Roche et al. 19331 ; [infante fc Pritchet 1995| [Brainerd, Small, k Mould 19351 ). 



Second, the typical mass of these objects (10 9 — 10 10 M Q ) is close to that of a dwarf galaxy 
[cf. , e.g., Babul & Rees (1992) scenario]. Finally, the number density of the clouds is 
reasonably close to what is observed for faint blue objects. 

Finally, we caution that it may be necessary to extend the simulation's dynamic range to 
1000 — 3000 cells in each dimension, in order to have both a larger box (L ~ 20 — 30/i _1 Mpc 
comoving) and a finer resolution (A/ ~ 10 — 20/i~ 1 kpc comoving) to ensure both a "fair" 
sample and fully resolved structures of the objects in question, while maintaining a sufficient 
mass resolution. 
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Fig. 1. — displays three dimensional isodensity surfaces at p cut = (3, 10, 30) (in panels a,b,c), 
respectively, at z — 3. 

Fig. 2. — Figure 2a shows two randomly selected slices of size 5 x 5/i~ 2 Mpc 2 with thickness of 
175/i -1 kpc at z = 3 (all lengths are in comoving units). The two left panels show the density 
contours with p cnt 10 < ^~ 1 - ) / 2 ,2 = 1,2,3..., and the two right panels show the identified clouds 
in the same slices, where two symbols are used to show two types of clouds: the filled dots 
represent "large" clouds with their sizes roughly indicating the actual sizes of the clouds, 
and the open circles are for the "small" clouds (see text for definitions). Figure 2b is the 
same as Figure 2a but with p cut = 10 for two different slices. 

Fig. 3.— shows the cumulative distributions of the cloud sizes at z = 3 for three different 
cases: p cut = (3, 10,30), weighted by A (thin curves) and by mass (thick curves). 

Fig. 4. — shows the redshift evolution of the median size of the clouds weighted by A (thin 
curves) and by mass (thick curves). Note that here the size (vertical axis) is in comoving 
length units. 

Fig. 5. — shows the distributions of the shapes of the clouds in the c/b-b/a plane for clouds 
at z = 3 at p cut = 3 (5a), p cut = 10 (5b) and p cut = 30 (5c). Each dot represents one cloud 
and the contours indicate the density of clouds (weighted by A) in the plane. The contour 
levels are incremented up linearly from outside to inside. We see that the most common 
cloud at p cut = 3 and 10 is an ellipsoid with axial ratio of ~ 1 : 2 : 4. At p = 30 there are 
two major concentrations of clouds at (b/a,c/b) = (0.35,0.85) and (0.75,0.85), indicating 
the existence of two distinct populations: filaments and near-spherical ellipsoids. 

Fig. 6. — shows the shape r\ as a function of size S at z = 3 for p cut = (3, 10, 30) (from top 
panel to bottom panel), respectively. There is a trend, with large scatters, that larger clouds 
tend to be less spherical and more filamentary. 
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Fig. 7. — represents the A-weighted cumulative distribution of rj for three cases at z — 3. 

Fig. 8. — shows the redshift evolution of the median shape paramter f] me d- 

Fig. 9.— shows the fraction of lines which are coincident at z — 3, N co>tot /N tot (where 
N co ,tot = N COtCommon + N COtdustering ), as a function of Ar. Six curves are shown for three 
values of p = (3, 10, 30) each at two values of D — (50, 150) km/s. 

Fig. 10. — shows the redshift evolution of N COyto t/Ntot at Ar = 100/i" 1 kpc in proper units, 
for four cases with (N HI ,D) of (3 x 10 13 , 150), (3 x 10 13 ,50), (3 x 10 14 , 150), (3 x 10 14 ,50). 

Fig. 11. — shows the ratio of coincident line pairs, each of which shares the same cloud, to 
the total number of coincident line pairs, N COtCommon / N COtto t, at z — 3. Six curves are shown 
for p = (3, 10, 30) with D = 50 km/s and 150 km/s. 

Fig. 12.— shows rms difference of coincident line pairs along the line of sight, -Rii rms , as a 
function of the sigthline separation perpendicular to the line of sight, Ar, for D = 50km/s 
and 150km/s at p cu t = (3,10,30). Also shown as the two long-dashed horizontal lines are 
the rms separations for the two D values, if clouds are small (compared to perpendicular 
separation) and randomly distributed. 

Fig. 13. — shows the redshift evolution of N COtCommon /N COttot for Nhi = 3 x 10 13 cm -2 (panel 
a) and Nhi = 3 x 10 14 cm" 2 (panel b), respectively. Four cases are shown in each figure 
with two separations, Ar = 50/i _1 kpc and Ar = lOO/i^kpc both in proper units, and with 
D = 50km/s and 150km/s. 

Fig. 14. — shows the correlation of Lya clouds at z = 3 for p cut = (10, 30). 

Fig. 15. — shows the redshift evolution of correlation length r (solid curve) defined as the 
length where correlation is unity (in comoving length units). Also shown as the dotted curve 
is the correlation length computed using the eight subboxes to illustrate the effect of the 
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simulation boxsize on the correlation. 
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Table 1. Cloud Statistics 



Redshift p cut # of Clouds Mass fraction (r sp ) (Mpc proper) (m c ioud) ( M Q ) 
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